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ABSTRACT 

We have carried out a search for the 2.14 pm spectroscopic signature of the close 
orbiting extrasolar giant planet, HD 179949b. High cadence time series spectra were 
obtained with the CRIRES spectrograph at VLT1 on two closely separated nights. 
Deconvolution yielded spectroscopic profiles with mean S/N ratios of several thousand, 
enabling the near infrared contrast ratios predicted for the HD 179949 system to be 
achieved. 

Recent models have predicted that the hottest planets may exhibit spectral sig- 
natures in emission due to the presence of TiO and VO which may be responsible for 
a temperature inversion high in the atmosphere. We have used our phase dependent 
orbital model and tomographic techniques to search for the planetary signature under 
the assumption of an absorption line dominated atmospheric spectrum, where T and 
V are depleted from the atmospheric model, and an emission line dominated spectrum, 
where TiO and VO are present. 

We do not detect a planet in either case, but the 2.120 pm- 2.174 urn wavelength 
region covered by our observations enables the deepest near infrared limits yet to be 
placed on the planet /star contrast ratio of any close orbiting extrasolar giant planet 
system. We are able to rule out the presence of an atmosphere dominated by absorption 
opacities in the case of HD 179949b at a contrast ratio of F p /F* ~ 1/3350, with 99 
per cent confidence. 

Key words: Line: profiles - Methods: data analysis - Techniques: spectroscopic - 
Stars: late-type - Stars: individual: HD 179949 - Stars: planetary systems 



1 INTRODUCTION 

Recent models (|Burrows et al.l 120081 ) have focused on the 
role of a parametrized stratospheric absorbing molecular 
species which leads to abundant molecules such as H2O ap- 
pearing in em i ssion in certain wavelength bands. Similarly, 
IFortnev et ail {2008) have more specifically investigated the 
effects of TiO and VO which absorb much of the incom- 
ing stell ar radiation hi g h in t he atmospheres of the hottest 
CEGPs. IFortnev et ail (120081 ) expect that HD 179949b be- 
longs to this class of hot CEGPs where the absorbing species 
lead to a temperature inversion and the formation of a 
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stratosphere. It is claimed that such models give better 
agreement with multi-wavelength mid-infrared Spitzer ob- 
servations for those CEGPs which ar e most highly i rradi- 
ated. Systems such as HD 209458b l|Burrows et al.l 120071 ; 
iKnutson etafl 120081 ) . HD 149026b and v And are among 
those systems which possess atmospheres most consistent 
with the presence of temperature inversions. 

In the near infrared, the strong 2.2 u.m bump due 
to the presence of strong H2O and CO molecular bands 
for 2.2 pm < A > 2.2 pm may not be presen t in s ystems 
which exhibit a stratosphere dBurrows et all [20081. This 
possibility was indicate d by observations of HD 209458b 
(Richards on et all I2003T ) which failed to detect the 2.2 pm 
bump in the atmosphere of HD 209458b. The weak H2O and 
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CO absorption transitions close to 2.2 [im are rather seen in 
emission while the 2.2 u.m F v /F* flux ratio is lower and 
natter owing to a spectral energy distribution which more 
closely resembles that of a blackbody. 

1.1 HD 179949 and its planet 

The presence of a close orbiting planetary system around 
the M v = 4 .09 ± 0.04 F8 dwarf star, HD 179949 w a s first 
reported bv lTinnev et all (j200ll ). lEggenberger et al.1 (120071 ) 
included HD 179949 in a search for stellar duplicity but 
found no companion star, indicating it to be a single star 
system. Further, a survey of stars known to harbour planets 
(Santos et al.1 12004! ) revealed that the [Fe/H] = 0.22 ± 0.05 
dex for HD 179949b places it in the most common metal- 
licity band, in a distribution significantly more metal rich 
than than for stars which are not known to harbour plan- 
ets. The reported velocity amplitude induced by the planet 
on its parent star, at K = 101.3 ± 3.0 km s _1 , with orbital 
period of 3.093 ± 0.001 d and orbital radius of 0.045 ± 0.004 
placed it am o ng th e closest orbiting planets. More recently, 
iButler et al.1 (|2006h have published revised system param- 
eters as a result of further monitoring of the system, find- 
ing K = 112.6 ± 1.8 km s~\ P = 3.092514 ± 0.000032 d 
and a = 0.0443 ± 0.0026 AU in a low eccentricity (e = 
0.022 ± 0.0 15) orbit. By studying emission in the Ca II H 
& K lines, IShkolnik et al.l (|2003l ) found evidence for per- 
turbation of the stellar magnetic field with a periodicity 
which coincides with the orbital period of HD 179949b while 
IWolf fc Harmanecl |2004) found the stellar rotation period 
of HD 179949 to be independent of the planetary orbital pe- 
riod, with Prot = 7.06549±0.00061 d. With no transit of the 
system reported, a minimum mass of Msini — 0.916 ±0.076 
Mq is found. The time of inferior conjunctio n is given by 
an ep h emeris of HJD = 2 4 51001 .510 ± 0.020 d jButler et al.l 
l2006l ). IWittenmver et al.1 l|2007l ) found no evidence for long 
period objects in the HD 179949 system and reported up- 
dated system param eters which are consistent with those of 
IButler et all (120061) at t he 1-cr level. 

ICowan et al.l l|2007h presented mid-infrared Spitzer ob- 
servations of the HD 179949 system which revealed a phase 
dependent light curve in phase with the planet's orbit and 
with a relative peak-to-trough amplitude of 0.00141 at 8 urn. 
This implies that less than 21 per cent of the incident stellar 
radiation is recirculated to the night side of the planet and 
contrasts with the other systems in their study, HD 209458 
and 51 Peg, which did not reveal photometric variation, sug- 
gesting a higher level of redistributi on of incident energy . 
While these observations were used bv lBurrows et al] ( 2008) 
to compare model fits, the data carry sufficiently high uncer- 
tainties that they enable a number of models with varying 
degrees of heat redistribution, absorber opacity and incli- 
nation to give reasonable fits. A hint that this degeneracy 
may be partially broken com es from the similarity b etween 
HD 179949b and v And b (jHarrington et al]l200d ). both 
showing only a small shift between the superior conjunc- 
tion ephemeris and lightcurve maximum. This may indicate 
that heat is re-radiated with out being carried downstream , 
as suggested by the results of lHarrington et all (|2006l . 120071 ) , 
due to the presence of a hot stratosphere which inhibits ad- 
vection in the lower atmosphere. 

In this paper, we present 2.14 \im observations taken 



with the Cryogenic high-resolution InfraRed Echelle Spec- 
trograph (CRIRES) at the Very Large Telescope and search 
for the signature of a planet which exhibits either absorption 
or emission features. We first describe observations and data 
extraction procedures (ij2| before describing the method (ij3| 
which involves searching for the faint planetary spectrum in 
a mean spectrum subtracted timeseries of spectra. A Gaus- 
sian matched filter mimics the radial velocity motion and 
varying strength of the planetary spectrum which is used to 
search for the best fitting model. The results are presented 
in ij4]and discussed in fj5] 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of HD 179949 were made with CRIRES at 
the Nasmyth focus of the Very Large Telescope, Unit 1 
(VLT1, Antu), on 2007 July 26/27 and August 03 and are 
presented in Table \T\ Four 1024x1024 InSb Aladdin-3 ar- 
rays were used to obtain each spectrum in the wavelength 
ranee AA = 21215 - 21740 A. Observations were made in 
ABBA nodding sequences to enable subtraction of sky back- 
ground and lines and to remove the effects of detector de- 
fects such as hot pixels. The recorded frames of HD 179949 
comprised 4 x 25 s exposures (see Table [T}. Observations of 
early type Standard stars were also made to enable moni- 
toring of changes in telluric line strength throughout each 
night. For the first and second nights, the seeing varied in the 
ranges 0.39" -0.97" and 0.40" - 1.4" respectively. The seeing 
was generally ~ 0.5" for observations of HD 179949, while 
the humidity was very low at 3-8 per cent and 12-20 per 
cent on each night respectively. 

2.1 Data extraction 

The worst cosmic ray events were removed at the pre- 
extraction stage using the Starlink FIGARO routine bclean 
l|Shortridgelll993l ). The Aladdin-3 arrays used with CRIRES 
exhibit an odd/even non-linear sensitivity pattern between 
adjacent columns (i.e. in the cross-dispersion direction) 
on dete ctors 1 and 4 and adjacent r ows on detectors 2 
and 3 (|Siebenmorgen fc Smettel 120081 ). The effect is not 
removed by simple flatfielding or nodding between A and 
B positions and failure to account for this effect results 
in extracted spectra which display a saw-toothed pattern 
on a pixel-t o-pixel scale. We implemente d the method 
described bv lSiebenmorgen fc Smettel {2008) to account for 
this effect before further extraction of the science frames 
was carried out. Each spectrum was then extracted by first 
subtracting the neighbouring frame in and ABBA sequence 
to eliminate sky background and detector cosmetic effects. 
This method means that during extraction of the spectra, 
it is not strictly necessary to model the sky background. 
However, we fitted low order polynomials (straight line) 
across the profile at each dispersion position to account for 
any slight changes in the sky background from exposure 
to exposure. Due to the dead space between detectors, 
it is not possible to record a continuous spectrum with 
CRIRES over a selected wavelength range. We therefore 
carried out the extraction process in the same manner as we 
would for four individual orders of a cross dispersed echelle, 
with the relevant gain and readout noise for each detector. 
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Table 1. CRIRES/VLT1 observations of HD 179949 for UT 2007 July 26 and August 3. 



Object 


UT start of 


UT start of 


Time per Number of 


Number of 


Comments 




first frame 


last frame 


exposure [sees] co-adds 


per frame 


observations 




UT 2007 July 26/27 


HD 179949 


23:51:46 


00:06:52 


25 


4 


8 


Main F8V Target 


HD 182645 


00:18:04 


00:24:35 


50 


2 


1 


B7IV standard 


HD 179949 


00:39:49 


03:26:40 


25 


4 


76 


Main F8V Target 


HD 158643 


00:40:07 


03:51:41 


10 


6 


4 


A0V Standard 


HD 179949 


04:04:28 


07:54:52 


25 


i 


104 


Main F8V Target 


HD 212581 


08:09:05 


08:35:48 


50 


4 


8 


B9.5V Standard 


HD 212301 


08:48:28 


09:51:06 


50 


4 


16 


F8V target 


HD 212581 


10:05:19 


10:10:31 


12 


6 


4 


B9.5 Standard 


UT 2007 Aug 03 


HD 158643 


00:18:55 


00:23:17 


15 


4 


4 


A0V Standard 


HD 179949 


01:08:19 


02:44:18 


25 


4 


41 


Main F8V Target 


HD 158643 


03:05:54 


03:17:00 


10 


6 


8 


A0V Standard 


HD 179949 


03:31:11 


07:22:23 


25 


i 


104 


Main F8V Target 


HD 212581 


07:39:34 


07:53:31 


15 


6 


8 


B9.5V Standard 


HD 212581 


07:56:14 


08:36:11 


20 


6 


16 


B9.5V Standard 


HD 212301 


08:49:44 


09:35:04 


50 


4 


12 


F8V target 


HD 212581 


09:44:38 


10:24:05 


20 


6 


16 


B9.5V Standard 



The typical gain for a detector is 7 electrons/ ADU while 
the r.m.s readout noise is 10 electrons. The spectra were 
extracted using echomop's imple mentation of the optimal 
extraction algorithm developed bv iHornel ljl986t ). echomop 
propagates error information based on photon statistics 
and readout noise throughout the extraction process. 



2.2 Wavelength calibration 

The wavelength range of our observations was chosen as a 
compromise between minimising the effects of tellurics, max- 
imising the expected number and strength of planetary spec- 
tral lines (See §4) and maximising the expected planet/star 
flux ratio. Calibration frames using ThAr lamps were taken 
during the observations, but only yielded typically one or 
two lines per order in the observed spectral ran ge. We there- 
fore fol lowed the successful approach used in iBarnes et al.l 
l|2007al lbl). where telluric features in the observed spectrum 
of the standard star were used. A spectrum generated from 
a HITRAN line list (|Rothman et al.i r2005) enabled identi- 
fication of observed features with theoretical wavelengths. 
A quadratic function was sufficient to describe the non- 
linearity between pixel position and the wavelength scale, 
with between 8 and 23 telluric features used for the spectrum 
on each of the four detectors. A residual r.m.s. of < 0.042 
times the spectrograph resolution element (a 0.4" slit is 
equivalent to 0.429 A at the mean wavelength of the ob- 
servations ) was obtained resulting in wavelength ranges of 
21214.9-21330.3 A, 21363.0 - 21473.9 A, 21503.5-21610.2 A 
and 21638.3-21740.0 A respectively being recorded on each 
detector. 



3 DATA ANALYSIS 

3.1 Removal of stellar and telluric spectra 

In order to maximise our chances of detecting a signal we 
first remove contributions from sources which contaminate 
the planetary spectrum. We have found that the most ef- 
fective means of removing the stellar and telluric lines in- 
volves subtraction o f a scaled and shifted master spectrum 
(IBarnes et al ] |2007bl ). For each night of observations, a mas- 
ter spectrum is created by cross-correlating and shifting all 
observed HD 179949 spectra to the nearest pixel. The re- 
sulting spectrum contains an average of the telluric, stellar 
and planetary signal with the planetary signal smeared out 
due to its radial velocity variation with orbital phase. The 
phase ranges over which we chose to observe HD 179949 en- 
sure that any planetary absorption signature will be Doppler 
shifted from -113 km s _1 to +91 km s _1 during the planet's 
orbit about HD 179949. 

After subtraction of the master spectrum, the remain- 
ing principal variation from spectrum to spectrum is seen as 
a changing telluric line strength throughout the night due 
to changes in humidity and airmass of the observed target. 
The telluric lines thus change their strength relative to the 
stellar lines, and relative to each other. The consequence 
of the latter effect is that accurate removal of time varying 
telluric features can only be carried out effectively when in - 
di vidual lines ar e resolv ed, as shown by Bai ley et al.l (|2007h . 
In IBarnes et al.l l|2007bl ). we fitted the master spectrum to 
each HD 189733 spectrum in turn by taking derivatives of 
the spectra and using splines to calculate the scale factor 
at points across the spectra. This process can account for 
those lines th at behave independent l y ove r the night and 
is described in lCollier Cameron et all l|2002h (Appendix A). 
Care must be taken not to use too few scaling points, which 
would not adequately account for relative changes in telluric 
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Figure 1. CRIRES wavelength coverage of HD 179949. The four segments (left to right and top to bottom) are the spectral regions 
recorded on detectors 1-4 respectively. Plotted (top to bottom) in each panel are the mean spectrum for 2007 July 26, model telluric 
spectrum, model planet spectrum with TiO/VO removed (model 1) and model planet spectrum with TiO/VO stratosphere (model 2). 
For model 2, the transitions are in emission and are weaker relative to the normalised continuum in the 2.120 um 2.174 \im spectral 
window. The model spectra have been convolved with a Gaussian to match the instrumental resolution. The tick marks represent opacity 
wavelengths (before convolution) for each planetary model. 



line strength. Conversely, use of too many scaling points 
could effectively yield a close or perfect match between the 
individual and master spectrum since the continuum in the 
master spectrum could be scaled to mimic a weak absorp- 
tion (e.g. planet) feature. This would result in attenuation or 
removal of a potential planetary signal. We therefore do not 
include more scaling points than twice the effective spectral 
resolution. With these limitations, the effectiveness of our 
spectrum removal procedure is dependent on spectral reso- 
lution , and at the resolution of R ~ 15, 000 |Barnes et al.l 
l2007bl ) significant residuals inevitably remained. In contrast, 
the present CRIRES data have R ~ 50, 000, which enable 
more accurate telluric subtraction. 

Some tellurics nevertheless still remain at the 2-er level 
(JU after subtraction of the scaled master spectrum. We 
therefore reduced the strength of remaining systematics with 
principal components analysis of the spectral timeseries us- 
ing th e method first described in ICollier Cameron et al. 
1 20021) (Ap pendix B) and discus sed in both llBarnes et^rt* 
(2007a]) and lBarnes et all (|2007bh . The data at fixed wave- 
length positions are decomposed into a series of eigenvectors 



which the main orthogonal trends in the timeseries of ob- 
servations. Removal of the first four principal components 
ensure that the remaining systematics due to time varying 
telluric residuals are subtracted while not significantly at- 
tenuating a moving planetary signal. At the high S/N ratios 
(see of the observed spectra, complete reduction of 

residuals to levels close to 1-a is not possible. 

3.2 Deconvolution and Atmospheric Model 

To enable absorption or emission features to be detected, 
we deconvolve our observed spectra using a model of 
the planetar y atmosphere . In e ffect, a least squares de- 
convolution l|Donati et al.l 119971 ) is applied to each ob- 
served spectrum by using the wavelength positions and line 
dep ths of all opacities in the wavelength range o f inter- 
est (|Collier Cameron et"ai"1ll999l ; lBarnes et al.ll2007al lbl). For 
HD 179949b, we model an atmosphere with solar metal- 
licity for a temperature, T = 1250 K and surface gravity, 
log g = 1.33 ms _1 . The effective wavelength of the single 
mean profile is determined by the distribution of line opac- 
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Figure 2. Orbital phase diagram for HD 179949b. Phase 0.00 
corresponds to the time of inferior conjunction. The large filled 
circle at centre represents HD 179949 while the small filled circles 
represent the phases of observations. 



and (2) where TiO and VO lead to formation of a strato- 
sphere, resulting in formation of weak lines in emission (here- 
after model 1 and model 2). The form of the spectrum in 
model 2 is much more blackbody-like than in model 1 and 
results in an overall lower value for the flux ratio in the 
wavelength region of our observations. For a fuller descrip- 
tion o f the model opacities and setup see Ferguson et al.l 
[|2005h . lBarman et al.l <[200ll ) and lBarman et al.l 120051 ). The 
re sulting models are ex pect ed to be very simila r to those 
of lBurrows et al.l (|2008h and lFortnev et all <[200Sl ) described 
above. Fig.[T]shows the mean normalised spe ctrum from July 
26, th e model normalised telluric spectrum l|Rothman et al.l 
2005) and the normalised planetary spectra for models 1 and 
2. 

After co-addition of ABBA sequences, a total of 46 spec- 
tra were observed on July 26 and 27 spectra were observed 
on August 2. The mean S/N ratio, measured from residual 
ABBA spectra after removal of the template and principal 
components, was 959 ± 57 and 893 ± 73 for the night of 26 
July and 02 August respectively with a combined S/N ra- 
tio of 935 ± 71. Deconvolution using model 1 yielded S/N 
= 8263 ± 721 while model 2 yielded 2589 ± 252, indicat- 
ing effective gains of 8.84 and 2.77 respectively. The lower 
S/N achieved through use of model 2 is a consequence of the 
weaker emission lines (see Fig. [TJ ■ 



ities, op acity strengths and number of counts at each wave- 
length [[Barnes et all ll99Sn . The mean weighted effective 
wavelengths of the deconvolved profiles are A c = 2.1434 u.m 
and A c = 2.1428 \im for model 1 and model 2 respectively. 
The similarity is a reflection of the uniformity in strength 
and distribution of opacities in the two scenarios. 

We have discussed the sensitivity of the model to 
mismat ches between in th e observed and model spec- 
tra in Barnes et alJ (|2007al lbh. For detection of absorp- 
tion/emission features it is necessary that a model with 
opacities which are as up to dat e as possible a re use d. Our 
models are generated using the iBarber et al.1 ([20061 ) water 
line list, H2O being the dominant species in the wavelength 
region of our observations. We have discussed the role of ab- 
sorbing species (J[J in the upper atmospheres of those plan- 
ets which absorb the highest amount of incident solar radia- 
tion. For cooler atmospheres, Ti and V are believed to rain 
out of the upper atmosphere, thereby giving no contribution 
to emergent spectra. In the case of the hottest CEGPs how- 
ever, T and V are present as oxides, namely TiO and VO, 
which absorb incident radiation high in the atmosphere lead- 
ing to heat ing and formati o n of a t emperature inversion . It is 
claimed bylBurrows et al.1 (|200St ). [Burrows et all (|2007n and 
iKnutson et all (|2008h (HP 209458 b) that multi-wavelength 
mid-infr ared observations are in better agreement with this 
scenario. iFortnev et al.l (|2008l ) have divided cooler and hot- 
ter planets into tw o regimes which they c all p L and pM class 
respec tively. Both iBurrows et alJ ([20081 ) and IFortnev et al.l 
(2008) discuss whether HD 179949 belongs to the former or 
latter of these classes, but given the relatively early spectral 
F8-F9V spectral type of HD 179949 itself, it seems likely 
that HD 179949 belongs to th e pM class. 

Building upon the work of lBarman et al.1 (|2005l ). model 
spectra are generated for the two cases (1) where TiO and 
VO are removed from the atmosphere, leading to dominance 
of strong bands of H2O absorption (in the region of interest) 



3.3 Matched filter 

We aim to isolate the phase dependent planetary absorp- 
tion signature from the total light of the HD 179949 sys- 
tem. Unlike the stellar spectrum, the planetary spectrum 
is Doppler shifted due to its orbit and is expected to ex- 
hibit a phase variation in st r ength due to irradiation of 
its inner face (|Barman et al.l 120051 ). The time-dependent 
changes in Doppler shift and brightness are modelled by 
using a Gaussian matched filter which has been descri bed in 
both iBarnes et al.l |2007al ) and iBarnes et al.1 (|2007bl ). The 
method is a adapted from t he reflected light studies whic h 
are detailed in Appendix D of lCollier Cameron et al.1 ([2002). 
By modelling the variation in Doppler shifted velocity and 
planet/star flux ratio (F p /F*), we are able, in the pres- 
ence of a strong detection to determine both the maximum 
planet/star flux ratio (eo = F p /F*) and the velocity ampli- 
tude of the planet, K v . Fig. [2]shows the phases at which ob- 
servations were made. Phases close to <f> = 0.5 are favoured 
because the hot inner face is presented to the observer, re- 
sulting in the highest planet/star flux ratios. 



3.4 Signal calibration and orbital inclination 

The matched filter is calibrated by means of a simulated fake 
planetary signal with known eo and K p . The fake signal is 
added to the extracted spectra before further analysis, which 
involves removal of the mean stellar spectrum and residual 
noise patterns. As HD 179949 is not a transiting system, 
the axial inclination of HD 179949 is unknown. However, 
since Rossiter-McLaughlin effect me asurements for a num- 
ber of system s , including HD 209458 (iBundv fc Marcvll2000l; 
Queloz et all bood ; IWinn et al.1 120051; IWittenmver et all 
20051) and HD 189733 (|Winn et al.ll2006l ). have shown that 
the orbital plane of the planet closely matches that of the 
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axial inclination of the star, we carry out Monte-Carlo sim- 
ulations to determine the most probable velocity amplitude, 
K p , and the orbital inclination, i, of HD 179949b. Using the 
empirical orbital parameters discussed in section Sjl] we find 
most probable values of K v = 115 km s" 1 and i = 47.6°. A 
matched filter 2D search as described above ( H3.3J1 for the 
best fitting parameters is scaled to ensure the the planet is 
recovered at the simulated value of eo for the fake planet. 

The significance of the result is assessed using boot- 
strap statistical procedures. The data are re-ordered in 
a way which randomises the observation phases while 
preserving the effects of any correlated systematic errors 
(|Collier Cameron et al.l [2002^ . Essentially, the order of the 
observations is randomised in a set of 3000 trials, thereby 
scrambling any true planetary signal while enabling the data 
to retain the ability to produce spurious detections through 
the chance alignment of systematic errors. The least squares 
estimate of logio eo and associated \ 2 as a function of K v 
enable us to plot 68.4, 95.4, 99.0 and 99.9 per cent bootstrap 
limits on the strength of the planetary signal. 



4 RESULTS 

We present the results of our search for the 2.143 [im high 
resolution spectroscopic opacity signatures of HD 179949b 
in Fig. [3] The timeseries deconvolved spectra spanning both 
nights of observations are shown in Fig. [3] for model 1 (top 
left) and model 2 (bottom left). In the former scenario, 
since the opacities are expected to appear in absorption, 
a planetary signature would be seen as a dark trail while 
in the latter scenario, a planetary signature would be seen 
as a light trail relative to the mean background. In each 
case, the dashed line represents the velocity path of a plan- 
etary signature at the most probable velocity amplitude, 
K p = 115 km s _1 , and does not represent a detection. 

The dark regions in the \ 2 plots of logio eo vs K v repre- 
sent local enhancements in y 2 for particular combinations of 
parameters. These features are all present within the 68.3, 
95.4 per cent confidence (1-2 a) intervals and are due to 
systematic features which remain in the residual spectra. It 
is possible to rule out a number of features where the false 
alarm probability would be high. An estimate of the upper 
limit to the velocity amplitude of K p (max) — 155.8 ± 9.1 
km s _1 is determined from the orbital period and semi- 
amplitude. Planets detected at velocities much in excess of 
this limit can thus be ruled out. Similarly, a planet with low 
orbital inclination would be expected to yield a much lower 
planet/star flux ratio owing to the heated inner face of the 
planet appearing foreshortened to the observer. 

For model 1, we are unable to detect a planetary sig- 
nal at the most probable velocity amplitude at a level of 
logio eo = —3.67 with 95.4 per cent confidence (2-cr). This 
equivalent to F p /F* ~ 1/4600. Similarly, for model 2, we 
do not detect a planetary signal with logio eo = —3.30 
(Fp/F, ~ 1/2000) with 95.4 per cent confidence (2-cr). The 
implications and significance of these results are discussed 
in |g] 

The epoch of inferior conjunction and orbital period 
(phase, 4> = 0.0) of HD 179949 carry uncertainties which 
may result in significant drift of the phasing by the time of 
observations. A drift of up to ± lhr 18 min or (f> = ±0.0175 



could have accumulated from the ephemeris up to the mid- 
epoch of the present observations. Since such a drift in ob- 
servations may be responsible for our inability to detect the 
planetary signal, we investigated shifting the time of infe- 
rior conjunction by ±15 mins intervals in this range. No 
clear candidate signature was however detected at any of 
the shifted ephemerides. 



5 DISCUSSION 

The significance of our analysis is presented in Fig. U which 
shows the model and observed planet/star flux ratios. Fig. 
|3] (left) indicates the 99.9, 99, 95.4 and 68.3 per cent upper 
detection limits for an atmosphere which contains no high 
altitude absorbing species which lead to a temperature in- 
version. In this model, only absorption lines are present and 
our analysis is able to reject this scenario with 99 per cent 
confidence at a level of logioeo = —3.53 (F p /F t ~ 1/3350) 
for a planet with the most probable velocity amplitude. The 
model 2 scenario, where TiO and VO are able to absorb in- 
coming radiation high in the atmosphere and form a temper- 
ature inversion, leads to a spectrum resembling a blackbody 
with weak emission lines superimposed (Fig. [4] right). Since 
the lines are weak, we are unable to reliably reject this model 
with our analysis since the 95.4 per cent (2-cr) upper limit 
of logw eo = -3.30 (F p /F, ~ 1/2000) at the most probable 
velocity amplitude is higher than the model contrast ratio of 
logioeo = —3.46 (F p /F* ~ 1/2900) at the mean wavelength 
of the observations. 

Without a clear detection of the signal, it is difficult to 
constrain the model predictions further since the flux ratio 
is dependent on the as yet unknown orbital inclination of 
HD 179949b. The prior distribution for K v indicates that 
the relative chance of detecting the planet with _Kp(max) 
(i.e. i = 90°) is ~ 9.7 per cent of a detection at K p . Simi- 
larly, an orbit with i.e. i = 30° carries only a 4.7 per cent 
relative chance. These extremes are represented in Fig. 2] 
by the grey spectra which have simply been scaled by the 
flux expected from the planet at i = 30° and i = 90°. Even 
with i — 30° , we can rule out the presence of a planet with a 
model 1 atmosphere with 99 per cent (marginal) confidence. 
Conversely, for model 2, were i = 90°, we would marginally 
detect a planet with 99 per cent confidence. There is no 
evidenc e for eclipse e v ents in HD 179949b however. 

In H arncs et al l (|2007al lbh. we investigated the effect 
of model-observation line strength mismatch. However, the 
transitions of the str ongest H2O lines w hich contribute to 
the model spectrum l|Barber et al.|[2006T ). and hence the de- 
convolved spectra, are known experimentally, with signifi- 
cant uncertainties on only the very weakest lines. If we can 
rule out wavelength and oscillator strength uncertainties in 
the input H2O transitions as sources of mismatch between 
model and observation, uncertainties in line strength may 
arise from the planetary temperature structure which may 
be incorrect due to the assumptions of hemispheric-wide hy- 
drostatic equilibrium. 

Until a clear detection of a planetary signature 
is made at high resolution, a quantitative empirical 
model-observation mismatch constraint will not be pos- 
sible .^Nevertheless, successful detection of H2O at opti- 
cal (|Barmanl [2007) wavelengths, and of H2O and CH4 
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Figure 3. Left: Phased timeseries of the deconvolved residual spectra of HD 179949. The dashed sinusoidal curve represents the expected 
motion of a planetary signal based on the most probable velocity amplitude, K p = 115 km s _1 . For plotting purposes, all deconvolved 
profiles have been normalised to the same noise level to optimise the visible information content of the plot. Right: Relative probability 
X 2 map of planet-star flux ratio logio eo vs K p . The first 4 principal components have been subtracted from the timeseries spectra. The 
greyscale in the probability maps represents the probability relative to the best fitting model in the range white for to black for 1. 
Numerous enhancements in \ result from low level systematic residuals (see main text) which can not effectively be removed. Plotted 
are the 68.3, 95.4, 99 and 99.9 per cent (bottom to top) confidence levels for a planetary signal detection. No planetary signal is detected 
at the most probable velocity amplitude, K p , which is represented by the vertical dashed line. 
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via transmission spe c troscopy at mid-in frared wavelengths 
dTinetti et aT]|2007bl; [Swain et alj|2008h . have relied upon 
models i Tinetti et al.l 2007al ) which contain r eliable molec- 
ular o pacities such as the H2O transitions of iBarber et al.l 
l|2006h . The mid-infrared observa tions of the spectral e nergy 
distribution of HP 189733b by ITinetti et ail (|2007bh and 
I Swain et al.l ()2008h rely on relatively broad spectral regions 
where the match between observation and model appears to 
be very close. While transmission spectroscopy is sensitive 
to the planetary atmospheric signature at the terminator, 
the same opacities are included in the model we have used 
for our d aysid e average spectrum, and have been used by 
iBarmanl |2008), to infer the presence of H2O and CO in the 
atmosphere of HD 189733b. 

Our inability to detect H2O (the only species in the 
wavelength interval of our observations) in the atmosphere 
of HD 179949b, can be explained by the expected model 2 
scenario and lack of required sensitivity in our data. On 
the other hand, f or HD 189733b, we rejected the expected 
model 1 scenario ijBarnes et al.ll2007bh using the same tech- 
nique as presented in this paper. This may be explained by 
a lower than expected contrast ratio resulting from efficient 
flux redistrib ution at the d eeper layers probed by shorter 
wavelengths (IBarmanl l2008h . This could also apply to our 
rejection of the model 1 scenario for HD 179949b. However, 
given the earlier spectral type of HD 179949 when com- 
pared with HD 189733, and subsequent higher irradiation 
of HD 179949b, as alrea dy discussed, the model 2 scenario 
seems more appropriate (jFortnev et al.ll2008l : Ib urrows et al.l 
2008). It is worth noting that the deeper atmospheric lay- 
ers probed by the shorter wavelength observations in the 
model 1 scenario may be blocked by the presence of the high 
altitude absorber(s) in the model 2 scenario. Hence a non- 
detection of HD 179949b with higher S/N ratio observations 
could not be explained in the same way as for HD 189733b. 
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Figure 4. Model and 2.14 |am planet/Star flux ratios for HD 
179949b. Shown are models for model 1, absorption (left) and 
model 2, TiO/VO stratosphere, emission spectrum (right). Some 
apparent emission features are due to stellar absorption lines and 
result from dividing the planet spectrum by the stellar spectrum. 
The black crosses indicate the mean model flux over the wave- 
length of observations for the most probable inclination. Upper 
limits corresponding to the levels in Fig. [3] (99.9, 99 and 95.4 per 
cent) are marked by horizontal bars (from top to bottom) for our 
analysis using each model. The 99.9 per cent observational up- 
per limit is coincident with the mean flux level for the model 1 
model. The width of the top horizontal bar in each case (99.9 per- 
cent confidence) represents the wavelength range of the CRIRES 
spectroscopic observations. For model 1, we can marginally ex- 
clude the model with 99.9 per cent confidence at the most proba- 
ble inclination of 47.6°. For model 2 we are unable to exclude the 
model since the 95.4 per cent confidence level (2-scr) is above the 
model flux ratio and artifacts are present between 2-sc and 1-sct 
(68.3 per cent confidence marked by the arrow). 



6 CONCLUSION 

We have applied a planetary atmospheric model with two 
physically different conditions in an attempt to detect the 
planetary signature of HD 179949b. We are able to reject the 
scenario in which Ti and V are depleted from the atmosphere 
of HD 179949b in agreement with the predictions of recent 
models. We have also assessed our ability to test the models 
in which TiO and VO are present in the upper atmosphere, 
resulting in a temperature inversion and lines appearing in 
emission. We are however unable to detect a planet using 
this model, or to reject the model since our data do not 
achieve the necessary S/N ratios required. The 99 per cent 
confidence limit is 2 times higher than the expected flux ratio 
at the most probable velocity amplitude of the planet and 
would thus necessitate a further six nights of observation 
before this level of sensitivity could be reached. 

In order to reliably test the most recent models which 
include a stratosphere and weak emission, cross dispersed 
spectrographs which offer greater wavelength coverage will 
afford higher S/N gain and thus improve our chances of de- 
tection. Equally, study of transiting systems reduce the de- 
generacy considerably in any analysis since the inclination is 
known and will be close, or equal to 90°. The greater max- 



imum brightness of a transiting CEGP system will increase 
our chances of detecting and characterising the planet. 
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